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A RECURRENCE FORMULA FOR SHEAR- LAC PROBLEMS 

By Paul Kuhn 



SUMMARY 



The analysis of the bending action in "box beams with 
appreciable shear deformation of the flanges becomes very 
difficult in the general case of variable cross section 
and loading. This paper presents a convenient method of 
solving the problem "by the familiar method of dividing the 
beam into a number of hays that can be assumed to have 
constant cross, section and loading. Application of for- 
merly derived shear-lag formulas leads to a general equa- 
tion closely analogous in form to the well-known three- 
moment equation. A numerical example and two comparisons 
between calculation and experimental results are included. 



INTRODUCTION 



In the sheet- stringer combinations typical of present- 
day aircraft construction, the distribution of the stress- 
es is materially influenced "07/ the shear deformation of the 
sheet. A basic theory for this "shear-lag" problem was 
discussed in reference 1 and the analysis of single- 
stringer structures, such as shown in figure 1, was in- 
cluded for simple cases. Reference 2 showed that the ad- 
ditional mathematical difficulties introduced "by multi- 
stringer structures, such as shown in figure 2, could he 
circumvented hy using "substitute single-stringer struc- 
tures." 

The analysis of even a single-stringer structure, 
however, presents difficulties in the general case, in 
which the cross section as well as the loading may vary 

„ alon? the length of the structure. In reference 3 were 

discussed various methods of dealing with this problem by 
means of numerical integration methods. These methods 

» furnish a comparatively simple soluti on , ~cons~i dering the 

complexity of. the problem, but they are inadequate for 
taking into account all the possible variations of parame- 
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terg. One result of this weakness Is that the solutions 
may be reliable only near the root of the "beam, a serious 
drawback in all of the cases where it is not certain that 
the critical section will he near the root. For such 
cases, mathematically more rigorous methods of- solution 
are desirable. A principle for such a method was given 
in reference 1 in a very "brief form; the present paper 
shows in detail the most convenient manner of setting up 
the equations and also shows comparisons of theory with 
test data. 

BASIC ASSUMPTIONS AND . SIMPLIFI OAT I ON S 



For the sake of completeness, the "basic assumptions 
and methods of procedure will "be summarized. 

The "beam is idealized "by assuming each stringer to- 
gether with a certain effective width of sheet to be con- 
centrated at its centroid. The resulting idealized 
stringers are assumed to carry only longitudinal stresses. 

The sheet is assumed to carry all the shear. When 
the sheet buckles, allowances must "be made for the reduced 
shear modulus. When reasonably heavy stringers are at- 
tached to the sheet "by two rows of- rivets, it is suffi- 
ciently accurate to consider only the "free" width of 
sheet as being subjected to shear deformation. 

The transverse stiffness E of the cover is assumed 
to be • infinite. y 

The shear stiffness G of shear webs in beams is as- 
sumed to be infinite. 

For a beam with a number of stringers, which is the 
usual case, the stringers are replaced by a single "sub- 
stitute" stringer (reference 2) 

The material on the side opposite to the one being 
analyzed is, in general, assumed to be concentrated into a 
compact flange attached directly to the shear web, con- 
verting the. closed box into an open box. For example, 
when the tension side of the beam is being analyzed, all 
the compression material is assumed to be concentrated 
into a flange attached to the shear web. The effective 
centroid for this flange may be estimated, with some al- 
lowance being made for shear-lag action, on the basis of 
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experience or of a first calculation. Typical cross sec- 
tions of the resulting substitute single-stringer "beams 
are shown in figures 3(a) and 3(b). As indicated "by the 
dotted lines, the boams are assumed to "be symmetrical about 
a longitudinal plane so that only one-half of the structure 
needs to be analyzed, 

The substitute "beam gives the stresses in the corner 
flange of the actual multi stringer "beam. The calculation 
of the normal stresses in the stringers and of the shear 
stresses in the sheet of the actual "beam is dealt with in 
reference 2. 

All calculations given in this paper deal only with 
single-stringer beams. The subscript B denoting sub- 
stitute (reference 2) is therefore omitted as unnecessary. 

More detailed discussions of these . assumptions and of 
the basic theory, as well as the symbols and the sign con- 
ventions used, will be found in references 1 and 2. 



GENERAL PRINCIPLES OP ANALYSIS 



Direct analytical solutions of the shear-lag problem 
have thus far been given only for simple cases. Refer- 
ences 1 and 2 give solutions for beams in which the cross 
section is constant and the running shear in the web 
Sw/hfl- is constant. 

Actual beams will usually have a variable cross sec- 
tion and loading. Numerical solutions for such cases are 
obtained by a familiar expedient: The beam is divided into 
a number of bays so that no appreciable error is committed 
by assuming that the cross section and the running shear 
in the web are constant within each bay. The known analyt- 
ical solutions can then be applied to each individual bay. 

The complete numerical solution requires the knowl- 
edge of (r + l) boundary conditions for a beam with r 
bays. Two of these conditions are furnished by the known 
conditions at the ends of the beam, the root, and the tip. 
The remaining conditions are furnished by the principle of 
elastic continuity, which requires that the elastic defor- 
mation at the outboard end of any bay must equal the elas- 
tic deformation at the inboard end 'of the next bay adjoin- 
ing it. 
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DERIVATION OF EQUATIONS FOE BEAMS WITH FLAT COVER 
aND ' CONSTANT- V/IDTH AND DEPTH 



The beam is divided into bays as shown in figure 4, 
and the bays and the stations ate riumbered as indicated. 
Figure 5 shows the two adjacent cells n and. n + 1 as 
free bodies on~whiGh the forces act. For' the derivation 
ofr-tire ^equations, it will be assumed that the upper side 
of the beam is being analyzed for shear-lag effects; the 
lower side is therefore assura-ed to be concentrated into 
a flange attached to the shear web. 

In figure 6, the flange and the longitudinal stringer 
forces Fy and. F L are separated into two groups of 
f 0 r c e s : 

(1) The group of forces shown in fi"gur~o ~6(a"), do- 

noted by the superscript P indicating that 
these 'forces are calculated by the ordinary 
bending theory, which assumes that piano 
cross sections remain plane and (tacitly) 
that the shear modulus is infinite. 

(2) The group of forces X showh~~in figure 6(b) 

representing the changes in stringer forces 
caused by shear deformation of the cover 
sheet-. 

Since the first group of forces is in equilibrium by it- 
self, the second group must also bo in . equi li brium by it- 
self; that is,' the force Xy acting'on the flange at any 
given station must be equal and opposite to_ the force . X L 
acting on the longitudinal at the: same station. This re- 
sult was anticipated in figuro 6(b) .by omitting the sub- 
scripts .F and L from X. 

The chief reason for separating the forces into the 
P-^roup and the X-group is economy of arithmetic. The 
statically indeterminate calculation of tho shear-lag of- 
foots furnishes corrections only to the stresses computed 
by the ordinary bending theory and can therefore be. carried 
through with a lesser decree of accuracy than would be re- 
quired if the total forces were chosen as statically indo- 
terminato unknowns. 

The groups of forces acting on each bay are associat- 
ed with shear forces in the cover shoet. The deformations 
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caused "by tlieae shear forces will "be calculated, for the 
inboard end of "bay n and for the outboard end of "bay 
n + 1, If the two deformations are equated in accordance 
with the principle of elastic continuity, a general ex- 
pression will be obtained that will relate the X-forces at 
three successive stations with "the properties of the two 
bays between these stations. 

Figure 7 shows the cover of bay n with the forces 
acting on it separated into groups. The P-gfoug calculat- 
ed by the ordinary bending theory is shown in figure 7(a). 
The shear stress accompanying this group causes the shear 
deformation ^ni^ a * 'the inboard end of ' the bay as indi- 
cated, the superscript P indicating the force gr_oup caus- 
ing the deformation and the sub-subscript i_ 'denoting the 
inboard end. Figure 7(b) shows the X n -group acting at the 

inboard end of the bay. The shear stresses associated 
with this group cause the shear deformation Y ni a , Fig- 
ure 7(c) shows the X n _ x -group acting at the outboard end 
of the bay; the shear stresses associated with this group 

cause the shear deformation a-i a t the inboard end. 

The total shear deformation at the inboard end of bay n 
is 

Y =Y P +Y Xq + Y X n-»i (l) 
n i n i n i n i . 

Similarly, the . shear deformation at the onboard end of 
bay n + 1 is 

where the sub-subscript o denotes the outboard end. The 

SO 

familiar formula T = ~r of the ordinary bending' theory 

?iT6S _ - - 



Y * = Y * = Y P =/- » r .^V (3) 



i o ~ &e b w tA T & e 



The numerical 
for each bay, 
factors. 



values of this deformation can be computed 
using average values for the individual 
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By application of formulas (S-8) of reference 1, and 
introducing new 'symbols p and q_. f or conveni ence , it 
is found that 

^n> - -I„ r ~ Q ~ tanh(KL")^ = ~** Pn . (4) 
^i 1 *" 1 = X n-x = Xn . ;i qn ■ (5-) 



x n+ 1 „ K n+ 1 



Y (n+-i) 0 = ~ x n+i — — ~7ZT,7TTT~~ = " Xn+1 gn+1 



*n+i & e„ +1 sinh(KL) n+1 



In these formulas, A^ = Ay + A L , and the parameter Z 
is defined "by 



■g- 2 = lei (J 

1 b \A F • aj. 



\.Atp • At/. 



Substituting the expressions (4) to (7) into (l) and 

(2) and equating Y n± to Y (n+i) Q Si've the general re- 

lat ion 

P P 

X n-i In - x n ( ?n + Pn+i> + x n+i *n+i = ~ Y n + Y n+i 



which, is the "basic recurrence formula for the shear-lag 
problem. The unknown forces -X are obtained- by solving a 
system similar to the familiar system of three-moment equa- 
tions; i.e., a system of- n equations, each equation ex- 
cept the last one and usually the first one involving 
three unknowns. 
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BOUNDARY COHDITIONS 

1 

The first equation of the system is 

X oll - X x ( Pl "+ p a ) + X^ - - Y x + Y a (10) 

If only transverse, loads are applied to the shear web 
(fig. 8(a)), 



x o 



= 0 



If a couple P h ff is introduced at the tip of the shear 
web (fig. 8(h)), then 

x 0 =P^ (ID 
The last equation of the system is 

X r _ x q r - X rPr = -Y r + Y r+1 (12) 

In this equation, is tile "equivalent" shear defor- 

mation of the foundation to which the beam is attached. 
This deformation ie defined by 

Vi = i ■ - < 13) 

where. 8 is the relative longitudinal displacement of the 
root fittings at the flange I" and at the longitudinal L 

The conditions at the root will fall under one of 
three classifications-: 

(1) Rigid foun dation .- Although a perfectly rigid 
foundation is physically impossible, the foundation may 
be stiff enough to be considered rigid, and the displace- 
ment 8 will be zero. -In -practical cases, the condition 
of zero displacement may be obtained by symmetry of the 
beam about the plane of station r. 

(2) F oun dation yielding elasticall y o r inelasti call y 
The case of a foundation that yields elastically or inelas 
tically occurs, for instance, when the stresses in a wing 
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are carried through the f uselag" e "by t-ransverse members; 
the elongation of these members under load permits rela- 
tive displacement a 6 of the root fittings on the win?. 
These displacements' may "be obtained by successive approxi- 
mation; in the case of elastic yielding, a carry-through 
member may be treated as an additional bay of the beam, 
adding one unknown and one equation to the system, 

(3) Longitudinal free at the root .- The. case of the 
longitudinal that is unconnected at the root occurs, for 
instance, at the tension side of a beam when the stringers 
butt against the center section but are unconnected to it 
(or when the connection is too flexible to be taken into 
account structurally). In this case, the last equation 
of the system disappears, and X r is determined by static 

considerations to be 

i- = — ±* (14) 
T h'jy Arp 



BEAM WITH CAMBERED COVER AND TAPER IN 
DEPTH AND WIDTH 



The basic formulas given in the preceding sections 
require some modifications when the cover is cambered, when 
the beam is tapered in depth,, in width, or in both. 

In a beam with cambered cover as shown in figure 3(b), 
the parameter K is defined by (reference 2) 




(15) 



In addition to the X— group, it is necessary, for reasons 
of static equilibrium, to introduce also a group of forces 

X— — acting as shown in figure .9. These forces do not ap- 
h W 

pear in the statioally indeterminate calculation because 
they do not cause shear defrrrmati on of the cover? they 
must be taken into account, however, when the at-resses in 
the flanges are computed. 
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In a "beam with, a' flat cover, s'aear-la? action in the 
flat-cover side of the "beam does not ' change .the total 
force on the opposite side of the "beam at any given sta- 
tion; for instance, shear-lag action on the tension side 
does not change the total force on the compression. side 
if the tension side has a flat cover. Shear-lag action on 
a cambered cover, however, does change -the total force - on 
the opposite side of the "beam "by the amount Xt^— because 

the 'shift of force between the longitudinal and the flange 
on the cambered side reduces the effective depth of the 
beam. 

In a beam with cambered covers on both sides (fig. 
10) and with cross sections symmetrical about both axes, 
it is obviously unnecessary to analyze one side of the 
beam at a time; the conversion of the closed box into an 
open box by combining all the material oh one side into a 
single flange attached to the web is therefore omitted. 
The parameter K is defined by _.. 



Sb< V A ¥ A L j 

and the additional force group acting on the upper and the 
Z. v;cr comer g c 
lower corner flanges is X 

In a beam with cambered cover, the basic equation (3) 
must be written in the more general form 

& e I* &d » 

Also, equations (ll) and (14) for the boundary condi- 
tions must be written in the mo.re general form 

X Q = Ph w ^ — . _. (Ha) 

and 

X r = .. _. (14a) 



In a beam tapered in depth, the shear S in equations 
(3) and (17) must be understood to be the difference between 
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the -external shear and the vertical components of the 
stringer forces as given by the familiar formula 

S = S E - |j tan 6 (18) 

where h is the effective depth of the beam and 8 is 
the inclination between the effective compression f-lange 
and the effective tension flange. For a beam with cum- 
bered cover, equation (18) may be written in the more gen- 
eral form 

S = S E - ^ tan 8 (19) 

where CI is the static moment of Ay and A^ about the 
centroidal axis. 

Serious theoretical difficulties arise when the beam 
has taper in wridth. The method of analysis developed thus 
far may be considered dependent upon two distinct basic 
theories; the engineering theory of bending and the sim- 
plified theory of shear deformation in a skin- stringer 
combination. Both theories break down when the beam has 
taper in width, so that the development of an entirely 
rational theory is impossible without utilizing or devel- 
oping considerably more refined basic theories. 

As a temporary solution for engineering purposes, a 
semirational modification of the procedure for beams with 
constant width may be used. The actual width b and the 
actual thickness t of the sheet are replaced by a ficti- 
tious width b f and a fictitious thickness t f . The fic- 
titious width at any given station is taken as 

"bf = * - *o < 20 > 

where b is the actual width at the station and b Q the 
actual width at the root; the introduction of this ficti- 
tious width converts the tapered beam into a beam of con- 
stant width equal to the width at the root of the tapered 
b earn . 

The fictitious thickness t f _ at any station may be 
so determined that the shear stiffnesB of the fictitious 
sheet will equal the shear stiffness of the actual sheet 
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at that station. Because the shear stiffness of the sheet 
is proportional to t/b, this method of procedure would" 
require that ' .~ 

t f = t . - (21) 
b 

The use of this fictitious thickness then compensates for 

the use of a fictitious - width; at any station along the 

span, the shear-lag parameter .KL of the fictitious "beam 

with constant width equals the parameter K of the actual 

beam with tapered width. The coefficients P n , q n , and *Y n , 

as defined "by equations (3) to (7), are changed because 
they involve ' t independently of b. 

Considerations of the limiting case of a beam tapered 
to a point tend to indicate that this method of compensa- 
tion is inadequate, and compari son^'with the test s of 
N.A.C.A. beam 4 described in reference"!! snows rather ^poor 
agreement. More satisfactory agreement between. thi s test 
and calculation is obtained when the fictitious thickness 
is taken as 

t f - * (if) 2 .. - iiz) 

The use of this value for the fictitious thickness changes 
not only the coefficients p, q, and "V but also the shear- 
lag, parameter K at each station. 

Physically, the use of equations (2Q_) and (2l) may be 
interpreted as taking into account only the direct effect 
of variable width on the shear stiffness and neglecting 
all the effects of transverse components of forces caused 
by the inclination of the flange toward the center string- 
er. The use of equation (22) represents an attempt to 
take into account the effect of this inclination of the 
flange. 

COMPARISON BETWEEN CALCULATION AND TEST RESULTS 



Only two tests are available for checking the method 
of analysis developed. One already mentioned is that of 
N.A.C.A. beam 4 described in reference 2. The second one 
is described in reference 4. The beam used in the second 
test had a constant cross section, but the load was applied 
in the form of. several concentrated loads. Figure 11 shows 
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the experimental and the calculated results for N.A.C.A. 
beam 4; figure 12 shows the results for the beam of refer- 
ence 4. As a matter of some interest, the u-solution of 
reference 3, as well as the X-solution of this paper, is 
shown in these figures. 

The agreement "between te~st and calculation is not all 
that could be desired for N.A.C.A. beam 4. It is possible 
that some of the discrepancy near the root was caused by 
defective workmanship in one flange, which may be the rea- 
son for the decided drop in st-rese of one flange near the 
root. Regardless of whether defective workmanship in the 
flange caused the discrepancy, an experimental check of 
the specific influence of taper in width appears desirable. 

The agreement for the beam of figure 12 is very much 
better, particularly in the critical region near the root. 
The X-solution predicts correctly, at least qualitatively, 
the reversal in stress near the tip; the u-solution does 
not. The u-solution depends on the total shear at the 
root and gives a very good approximation near the root; it 
cannot take into account the manner in which the load is 
distributed along the span and, concomitantly, it cannot 
be expected to give results in very close agreement with 
the facts outboard of the root region, particularly if 
there are discontinuities of loading or cross section. 
Whether such disagreement is of practical importance de- 
pends, of course, on the individual circumstances. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley -Field, Ta. , November 16, 1939. 
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APP13DIX 
Numerical Example 



As a numerical example for .the' method of_ analysis de- 
veloped, the analysis of N.A.O.A. beam 4 will he given. 
The dimensions of this "beam are given in reference 2. The 
single-stringer "beam substituted for the actual "beam is 
calculated by the method described in reference 2 as fol- 
lows* - 

The compression flanges,.. the tension flanges r _ and the 
longitudinals are assumed to be concentrated at their re- 
spective centroids. The cover sheet, which is 0<0114 inch 
thick, is assumed to be fully effective in aiding the lon- 
gitudinals, and the adjacent strips of sheet are added to 
the longitudinals or flanges. The shear web is also as- 
sumed to be fully effective in bending and is replaced by 
concentrated flanges of cros s— sectional area l/& A w . The 

area of the substitute longitudinal is calculated by for~ 
mula (4) of reference 2: 

sinh E_b 

The substitute camber is taken as - - 




The subscript S is dropped as unnecessary for the present 
paper. With the geometrical properties of the single- 
stringer substitute beam thus defined and listed in rows 1 
to 5 of table I, the average section modulus l/zj, can be 
computed for each' bay and is listed in row 6. 

The effect of taper in width is allowed for by intro- 
ducing the fictitious width bf and the fictitious sheet 
thickness t f5 according to formulas (20) and (22); the 
shear-lag parameter is then calculated by the formula 
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and listed in row 7; KL is then "listed for each hay, L 
"being 19 inches for each "bay, and. the values of EL, sinh 
KL , and tanh KL are found in rows 8, 9, and 10, re- 
spe.ct i vely . 

The next step is the computation of the coefficients 
p, (j, and V for the system of X-equations. Some slight 
modifications were made' in this particular case to reduco 
the amount of arithmetic. The shear modulus & appears 
as a factor in all Coefficients and consequently can "be 



canceled. The thickness t f for all hays involves the 

thickness t as a common factor, so that t can he can- 
celed. With these modifications, the coefficients "become 



K 



tanh KL h Q ' 



q = 



sinh KL h 0 ' 



and are listed in rows 11 and 12. 

Before the coefficient Y is computed, the shear for 
each "bay must "be computed "by formula (18) or (19); tho load 
is a concentrated load of 250 pounds applied at the tip, 
and the resulting shear S is listed in row 13. With 
this shear S and by the use of the modifications previ- 
ously mentioned, the coofficient is 



and is listed in'' row 14. 

With these, coefficients , the Bystem of- oquations 
formed according to the ?enaral equation (9) is 



0.01354 Xi 
0.02068 X 2 
0.02908 X 3 
0.0383 X 4 



-0.0791 X x 
-0.0921 X 2 
-0.1070 X 3 
-0.1248 X 4 
-0.0373 X 5 



+0.01354 X s 

+0.02068 Xa 

+0.02908 X 4 

+0.0383 X 6 



0.04 
0.06 
0.02 
0.02 
•14.90 
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The solution of this system of equations ?ives the 
values of X listed in row 1 of table II. The flange 
stresses caused "by the X-group are 

x a f v h W / 

and are listed in row 2 of tattle II. This tabulation com- 
pletes the X-solution of the substitute beam. The total 
flange stresses cr in the actual beam are found by adding 
the stresses a"y_ to the flange stresses . Op, calculated 

for the actual beam by the ordinary bending theory, as 
shown in rows "3 and 4 of table II. _ Figure 11 shows the re- 
sults graphically. Strictly speaking, the stresses within 
each bay sho-uld be calculated by using formulas (B-8) of 
reference 1 or (A-f5 ) of reference 2; this refinement, 
however, is -seldom necessary. 
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TABLE I 

Basic Properties of Beam and Coefficients 
for Becurrenoe Formula 



Bay 


1 


2 


3 


4 


5 


1 


Aj 


-0. 1465 


0.1535 


0.1600 


0.1665 


0.1735 


2 


A L 


.330 


.395 


.460 


.525 


.590 


3 


■b 1 


6.71- 


7.93 


9.15 


10.37 


11.59 


4 


h w 


3. 00 


3.60 


4.20 


4.80 


5.40 


5 


c 


.615 


.725 


* .840 


.955 


1.070 


6 




1.523 


2.117 


2. 805 


3.59 


4.47 


7 


K 2 


.01392 


.00916 


.00643 


.00471 


.00357 


8 


KL 


2.240 


1.815 


1.520 


■1.304 


1.16-3 


9 


si.niL KL 


4.643 


2.988 


2.177 


1.7-06 


1.443 


10 


tanh KL 


.978 


.940 


.909 


.862 ■ 


.822 


11 


P 


.0365 


.0426 


.0495 


.05 75 


.0673 


12 


1 


.00770 


.01354 


.02068 


' .02908 


. 0383 


13 


S 


225.0 


187.5 


16.1.0 


140. 7 


125 .0 


14 


7 P 


14.76 


14.80 


14.86 


14.88 


14.90 
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TABLE II 
Flange lorcea and Stresses 



Stati on 


1 


2 


• 3 


4 


5 


1 X 


0.4 


5.0 


25.1 


89.4 


272.5 


2 a x 


3 


38 


185 


631 


1847 


3 Op 


1795 


2690 


2920 


3040 


3015 


4 o- 


1798 


2798 


3105 


3671 


4862 



^chnical Note Mo . 





N.X.C.A. Technical TS 0 te No. 739 



<s — 


-Is * 


ir - 






i 


— o- 


f— : 
1 


r 







J 

h 

6-1 



(a) 




00 



figure 3.- Typical cross sections of substitute single-stringer beams* 
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Figs. 4 





Figure 5.- Forces acting on individual bays. 
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(a) P - forces. 




Figure 6.- Division of forces into P - forces and X - forces. 
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Figs. 7,8,9,10 
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Figure 7.- Forces and deformations on cover of bay n. All forces and 
deformations are shown as positive. 





(a) (b) 
Figure 8.- Types of loading a.t tip bay. 
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Figure 9.~ Group of X - forces on 
cambered beam. 



Figure 10,- Beam with both 
sides cambered. 
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Fig. 11 
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Figure 11.- Flange stresses in tapered beam with, tip load. Experimental 
data from reference 2. 
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